Coordinate Representation of the One-Spinon One-Holon Wavefunction and 

Spinon-Holon Interaction 
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By deriving and studying the coordinate representation for the one-spinon one-holon wavefunction 
we show that spinons and holons in the supersymmetric t — J model with 1/r 2 interaction attract 
each other. The interaction causes a probability enhancement in the one-spinon one-holon wave- 
function at short separation between the particles. We express the hole spectral function for a finite 
lattice in terms of the probability enhancement, given by the one-spinon one-holon wavefunction at 
zero separation. In the thermodynamic limit, the spinon-holon attraction turns into the square-root 
divergence in the hole spectral function. 

PACS numbers: 75.10.Jm, 71.27. +a, 05.30.Pr 



I. INTRODUCTION 

Landau's Fermi liquid theory applies to interacting 
electron systems that can be adiabatically deformed to 
a Fermi gas. If the interaction is smoothly switched off, 
the spectrum of a Fermi liquid reduces to the spectrum 
of a noninteracting fermionic system. The excitations of 
a Fermi liquid are given by quasiparticles and quasiholcs. 
Although their lifetime may be short, it always becomes 
infinite at the Fermi surface pL As one is concerned 
only with energies close to the Fermi surface, Fermi liq- 
uid's picture applies to a wide class of correlated systems. 
Experimentally, Landau quasiparticles are observed as a 
resonant peak at the Fermi surface in the spectral density 
of states measured at fixed momentum. 

Nevertheless, there are several low-dimensional 
strongly correlated systems where Landau's picture 
breaks down. In Luttinger liquids the spin and charge 
degrees of freedom "separate" and the quasiparticles and 
quasiholes are no longer elementary excitations [^|J^] . The 
same phenomenon was discovered by means of Bethe- 
ansatz-like techniques in exactly solvable models, like the 
supersymmetric t — J-model with l/r 2 -interaction, which 
we study in this paper Q. 

Physically a particle or a hole injected in a strongly 
correlated chain breaks up into particles carrying spin, 
but no charge -spinons- and charge, but no spin - holons. 
The quantum numbers of particles and holes "fractional- 
ize" . Spinons and holons are the true elementary excita- 
tions of the strongly correlated chains j^^j] . When widely 
separated one from each other, spinons and holons prop- 
agate independently, in general with different velocities. 
Such a phenomenon is usually referred to as "spin-charge 
separation" . As a consequence of spin-charge separation, 
the quasiparticle peak disappears and is substituted by 



a broad spectrum. This was experimentally detected in 
ARPES experiments on quasi-lD samples (l]]. 

Using Bethe-Ansatz solutions of ID models, one can 
work out the energy of a many-spinon many-holon state. 
In the thermodynamic limit the total energy is the sum 
of the energies of each isolated particle. However, this 
does not imply that spinons and holons do not interact. 
In recent papers j^J^] we carefully studied spinon inter- 
action in an exact solution of the Haldane-Shastry model 
(HSM) a model for a strongly correlated ID sys- 

tem with no charge degrees of freedom. We showed that 
spinons interact, although in the thermodynamic limit 
the energy of a many-spinon solution is additive. We 
also showed that their interaction is responsible for the 
lack of integrity of spin waves against decay into spinons, 
the true low-energy excitations of the system. 

In this paper we generalize the formalism introduced 
in Ref. || to study the interaction between spinons and 
holons in an exact closed-form solution of the supersym- 
metric extension of the HSM: the supersymmetric t — J- 
model with 1 / r 2 -interaction (Kuramoto-Yokoyama (KY) 
model fjJI.) The KY- model is a system of electrons lo- 
cated at the sites of a circular lattice, where double occu- 
pancy of a site is forbidden by strong Coulomb repulsion. 
Charge hopping, Coulomb interaction and spin-spin anti- 
ferromagnetic interaction are all inversely proportional to 
the square of the chord between the corresponding sites. 
Charge vacancies at some sites (holes) may be created by 
filling the system with fewer electrons than the number 
of sites. In the KYM one is concerned with both spin 
and charge degrees of freedom. 

We investigate the basic features of the KY-model by 
employing a formalism based on analytic variables on the 
unit radius circle. As in the case of the HSM, it easier to 
construct and visualize the spinon and holon excitation 
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by using a real-space formalism than by using the Bethe- 
Ansatz formalism. Within our formalism we derive a 
"real space" representation of the one-spinon one-holon 
wavefunction as a solution of an appropriate equation of 
motion. 

By taking the thermodynamic limit of the exact solu- 
tion of the equation of motion, we find that the probabil- 
ity of finding a spinon and a holon at large separations 
from each other is independent of the distance between 
the particles, while it is greatly enhanced when the two 
particles are on top of each other, phenomenon we refer 
to as "short-distance spinon-holon probability enhance- 
ment" . In the thermodynamic limit, the spinon-holon in- 
teraction assumes the same form as the spinon-spinon in- 
teraction derived in [|| , although the corresponding equa- 
tions of motion are completely different. The physical 
interpretation is the same as in the case of the spinon- 
spinon interaction: a spinon and a holon do not interact 
when they are widely separated from each other, while 
they exhibit a short-range attraction at short separations 

To show the effects of spinon-holon interaction on 
the hole spectral density, Ah(uj,q), we exactly calcu- 
late the contribution to Ah(w, q) from one-spinon one- 
holon states, A^ p ' ho (o;, q). The spinon-holon interaction 
has important consequences on the functional form of 
A^ p ' ho ((jj, q). The probability enhancement causes the 
overlap between the wavefunction for the localized hole 
and that for a spinon-holon pair to be significant, al- 
though not enough to form a spinon-holon bound state. 
The corresponding matrix element is enhanced - despite 
the fact that the density of states is uniform at low en- 
ergy - so as to make the hole excitation fully unstable 
to decay into a spinon-holon pair. Taking the thermo- 
dynamic limit of our result, we show that, as the size 
of the system increases, spinon-holon interaction turns 
into a square root singularity at the one-spinon one-holon 
creation threshold. Correspondingly, A^' ho (uj,q) shows 
no Landau quasiparticle's peak, but it rather exhibits a 
sharp singular threshold, followed by a broad branch cut 

& 

The paper is organized as follows: In section II we 
shortly review the KY Hamiltonian and its supersym- 
metry; In Section III we introduce the ground state of 
the KY model at half filling and its representation as a 
function of analytic variables on the unit circle. At half- 
filling, the ground state is the same as the ground state 
of the HS-model - a disordered spin singlet. We will 
briefly review some properties of the ground state, al- 
ready discussed at length in Q ; In section IV we analyze 
the one-spinon solution and review its relevant proper- 
ties; In section V we focus on the one-holon solution and 
derive its relevant properties. In section VI we derive the 
action of TLky on the one-spinon one-holon states, the 
energy eigenvalues, the corresponding eigenvectors, and 



their norm; in Section VII we write the Schrodinger equa- 
tion for the one-spinon one-holon wavefunction, whose 
solutions are simple polynomials. From the behavior 
of the one-spinon one-holon wavefunction, we infer the 
nature of the interaction between spinons and holons: 
a short-range attraction. The physical consequences of 
such an interaction are discussed at length in Section 
VIII, where we rederive an exact expression for the con- 
tribution of one-spinon one-holon states to the hole spec- 
tral function in terms of the spinon-holon wavefunctions 
and rigorously prove that this contribution is completely 
determined by the spinon-holon interaction. In the ther- 
modynamic limit spinon-holon interaction turns into the 
square root divergence in the hole spectral function, ob- 
tained by jfl. In Section IX we provide our main con- 
clusions. 

II. KURAMOTO-YOKOYAMA HAMILTONIAN 

The Kuramoto-Yokoyama model is defined on a lat- 
tice with periodic boundary conditions. Sites are 
parametrized by the Y-th roots of unity, z a (a — 
1, . . . , N.) Strong electron repulsion forbids double oc- 
cupancy at each site. Therefore, sites can be occupied 
by an f or a |-electron, or they can be empty. The 
number of empty sites can be tuned by means of an ex- 
ternal chemical potential which fixes the total charge of 
the system. The Kuramoto-Yokoyama Hamiltonian |l2| ] 
is a generalization of the Haldane-Shastry Hamiltonian 
0,0, where also charge dynamics is taken into account. 
It takes the form: 

- \ ^{ciaCp*) + + np) - ^n a n p - ^ \p , (1) 

where the Gutzwiller's projector 

P = iT( 1 - c aT c il c «l c «T) - (2) 

a 

accounts for the no-double occupancy constraint. 
Site occupation and spin operators are given by 

n a — C^jCaf + C^|C a | (3) 

^ = ^E c -C' c -' ■ ( 4 ) 

cr,(T 

where r a , a = x,y,z, are Pauli matrices. An empty 
site corresponds to a charge- 1 hole that can tunnel to 
nearby sites by means of the same inverse-square matrix 
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element characterizing the spin exchange and the charge- 
repulsion term in TLky- 

Usual bosonic symmetries of a t — J-like model are 
total spin, corresponding to the operator S — J^a^a, 
and total charge, corresponding to the operator N — 
^2 a n a- The equivalence of energy scales for magnetism, 
charge transport and charge interaction, causes the KY 
Hamiltonian to be supersymmetric, in the sense that it 
commutes with the electron or hole injection operators 

As in the Haldane-Shastry hamiltonian, since the com- 
plex variable z lays on the unit circle (z* — z^ 1 ), the in- 
teraction is an analytic function of the coordinates, that 
is: 

\z a -zp\ 2 (z a ~zp) 2 

Throughout the paper we use the representation in terms 
of the analytic variables z a jl5|J|. This turns out to be 
very useful for describing the properties of spinons and 
holons in real space. 

III. GROUND STATE 

In this section and in the following one we discuss the 
ground-state and the one-spinon eigenstates of the KY- 
model at half filling. At half-filling the KY- model reduces 
to the Haldane Shastry (HS) Hamiltonian 

Both the ground state and the one-spinon eigenstates 
are the same as for the HS-model. Since in || we have 
already applied our formalism to study basic properties 
of the HS-model, here we will only briefly review the main 
results in view of their extension to states where holon 
excitations are present. 

A. Ground State wavefunction 

Let N be even. We first give the representation of the 
ground state |\Pgs) m terms of the z-coordinates and 
then derive its energy. |\&gs) is defined in terms of its 
projection onto the set of states with M = N/2 spins up 
and the remaining spins down. If z\, . . . , zm are the coor- 
dinates of the up spins, one defines the state \z±, ■ ■ ■ , Zm) 
as: \z\, ... , z M ) = njii S i Ua=i c lj°) where |0) is the 
empty state. The projections are given by: 

M M 

I'&sOi, ■ • ■ ,z M ) = X\{z - z k ) 2 \\ Zj ■ (6) 

j<k j=l 



^GS is a polynomial in the analytic variables z\, . . . , Zm- 
Its norm was first computed by Wilson jl6| by using the 
following identity: 

M 

c M = yi n>-*ii 4 

zi,...,zm i<3 




(7) 

Basic properties of ^gs follow in this section, together 
with their derivation. 

B. Singlet State 

The ground state is a spin singlet. I'J'gs) is annihilated 
by both S z and S~. S z \^ G s) = because |*gs) has an 
equal number of | and | spins, while 

N 

[S~-$ gs ](z 2) ...,zm) = ^2( z 2, ■ ■ ■ ,z M \S~ |*gs) 

a=l 



1=1 K a = l J 1 

since E«=i z i = N S to (mod N) . 

As ^gs is a spin singlet, it takes exactly the same 
form if expressed either in terms of the |-spin coordinates 
z±, . . . , zm, or of the |-spin coordinates, 771 , ... , T)m, that 
is: 

%gs(zi,...,z m ) = ^GsiVlf-iVM) ■ (9) 

Eq.(Q) is proved in Appendix A, where we derive the for- 
mulas to relate the representation of the states of the 
system in terms of |-spin coordinates to the representa- 
tion in terms of |-spin coordinates. 

In the thermodynamic limit, half-odd spin chains ex- 
hibit a gapless spectrum, although they are not al- 
lowed to order |i~7| . Accordingly, |\&gs) is a disordered 
spin liquid state, and the spin-spin correlation function, 
X(z a ) = (^gs\S^S~\^gs}/(^gs\^gs), Mis off with 
the distance as (— l) x /x, thus showing absence of spin 
order §. 

C. Ground State Energy 

At filling-1/2, I^gs) is the ground state of Hky, with 
eigenvalue: 
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TLky\^gs) — TLhsY&gs) 



(10) 



Eq.([10) has been derived originally by Haldanc and 
ShastryyQ. In @, we re-derived it by means of our 
own technique, consisting in substituting sums over spins 
on the lattice with derivative operators acting on the ana- 
lytic extension of ^gs{zi, ■ ■ ■ , Zm)- The Zj's are allowed 
to take any value on the unit circle. After computing 
the derivatives, we constrain them again to lattice sites. 
In this subsection we review our technique, in view of 
its generalization to the case where the filling is 7^ 1/2 
and the dynamics of the system is described by the full 
KY-Hamiltonian. 

Since [S^S^^gs]( z i, ■ ■ ■ , z m) is identically zero unless 
one of the arguments Z\ , . . . , zm equals z a , we have 

N S + S~ ~\ 

E ~ - f.12 \*GS](Z1, ...,Z M ) 



\z a - Zp | : 



M N 



E E I - 3 i2*C?s(zi, ■ • ■ , Zj-l, Zp, Zj+l, Z M ) 



zp 



N-2 M £+1 at f. Tl , \ 



(11) 



The coefficients Ai are calculated in ||. They are zero 
for N > I > 2. Therefore, Eq.diTI) can be rewritten as: 



M 

E 

3=1 



(N-l)(N-5) N-3 2 d 
Zj 2~ Z * dz~ 
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1 3 92 \ j ^Gs(zi,...,Z M ) 



iV(iV-l)(JV-5) 2zj 



24 



Zj Zk 



2z) 



M 7 2 



M 

- y — 



X* GS (zi, . . . ,Z M ) 



N 

T 



zm) ■ (12) 



The "Ising spin term" , on the other hand, provides: 



N 



iV(iV 2 - 1) 

48 



M . 

£zl*i-*kl a J 



(13) 



Adding up Eq.(|12j) and Eq.flL3|) provides Eq.(|10|). 

Using the factorization property of Ti-as, Shastry 
proved that [\Egs) is the actual ground state of TLky at 
1/2-filling p8[ . The same proof can be rephrased within 
our formalism, as discussed in Ref. ||. 

The crystal momentum of the state, q, is defined (mod 
2ir) by the equation: 

y GS (z 1 z 7 ...,z M z) = e iq -$ GS {z 1 ,...,ZM) , (14) 

where z=exp(2-7ri/iV) . From Eq.(|i"4|) q can be either 
or 7r, according to whether N is divisible by four or not. 
In the former case V^gs equals itself, when translated by 
one lattice constant, while it equals minus itself in the 
latter case. 

Other relevant properties of |^gs) are discussed in de- 
tail in H and will not be analyzed here. 



IV. ONE-SPINON WAVEFUNCTION. 

The elementary excitations above the ground state of 
a correlated ID electron system are not Landau's quasi- 
particles, but rather spinons and holons. Spinons have 
been identified as the elementary excitations of a spin- 
1/2 ID antiferromagnet. They can be thought of as lo- 
calized spin defects carrying total spin-1/2, embedded in 
an otherwise featureless disordered singlet spin sea. As 
spinon excitations appear in the KY-chain at any filling, 
we study spinons at filling- 1/2, when the KY- model re- 
duces back to the HS-model. One-spinon states appear 
as states of the chain with an odd number of sites. In 
this case, the minimum possible value for the total spin is 
1/2, and the state for a localized spinon at s is created by 
constraining the spin at s to be J., in a surrounding spin 
singlet sea |l(],[ll],|9| . Since correlations in the ground 
state are short-ranged, in the thermodynamic limit it 
makes no difference whether one begins with an odd or 
an even number of sites. Therefore, in the thermody- 
namic limit, there is no way to distinguish between chains 
with odd number of sites or chains with even number of 
sites. States with half-odd spin are alleged eigenstates of 
TLky at half-filling with an odd number of spinons and no 
holons. In this section we briefly review the one-spinon 
wavefunction and discuss its properties. This was first 
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studied by Haldane and Shastry [ pL 0|Jl 1[ | , and discussed at 
length in within the framework of the formalism of 
analytic variables. 



A. One-spinon spin doublet 

Let N be odd and M = (N — l)/2. The wavefunction 
for a localized spinon at s takes Haldanc's form [fK)|fll[: 



M 



M 



9f( Zl , z M ) = Uizj - s)z 3 Y[( Zl - z d f , (15) 

where z\ , . . . , zm denote again the position of the |-spins 
and s is the coordinate of a lattice site where the spin is 
fixed to be J,. 

By definition, \&| p is an eigenstate of S z with eigen- 
value -1/2. In order to prove that it is a spin-1/2 state, 
we need to show that S~ annihilates it. Indeed, per 
Eq.(||) we have: 



JV 



S7Vf = , 



(16) 



which proves that <3>| p is the |-spin component of a spin 
doublet. 

^f(zi, . . . , zm ) is a polynomial of degree less than 
N + 1 in each variables Zj. Therefore, we may again 
apply Taylor's expansion technique used to calculate the 
ground state energy. Doing so, we find: 



+ M(4M 2 -l)-^M 2 l>* 
6 2 



sp 

s ' 



(17) 



provided that A satisfies the following eigenvalue equation 
for $f = \[f{z j -s): 



M(M — 1) - a' 



ds 2 



N- 3 



M - s- 



d 



A$ sp 



(18) 



One- J, -spinon energy eigenstates are given by propagat- 
ing one-spinon plane waves: 

. . . , Z M ) = 1 X)(»*r*f ■ ■ ■ - (19) 



J.27T.O .iV-1 

2 ( iV } 



m) 



(20) 



with < to < (N - l)/2 and A = m((iV - l)/2 - to). 

As the total crystal momentum of the state |^^) is 
given by 



7f, P = Zn 



2tt 



(to - 



(mod 2-7r) 



(21) 



the total energy may be rewritten as: 
Hky\*Z) = {-J(^)(N + | - ^) + 



(22) 



that is, the sum of a ground-state contribution, plus the 
kinetic energy of the propagating spinon. The one-spinon 
dispersion relation is correspondingly provided by: 



£(0 



(mod 7r) , 



(23) 



As extensively discussed in |T^JTl|,^| , the one-spinon dis- 
persion relation shows the typical features of spinon ex- 
citations. It spans only the inner or outer half of the 
Brillouin zone, depending on whether N — 1 is divisible 
by 4 or not, which corresponds to the absence of nega- 
tive energy states, i.e., to the absence of "antispinons" . 
The spinon dispersion at low energies is linear in q with 
a velocity 



spinon — ^ u 



(24) 



The half-band of single elementary excitations for odd N 
are the only 5=1/2 states without extra degeneracies. 
The ground state of the odd-N spin chain is 4-fold de- 
generate and is given by l 1 ?^) for m = and (N — l)/2 
and their j counterparts. This corresponds physically to 
a "left-over" spinon with momentum ±7r. 

The spin density in the state ^> s s p as a function of the 
spinon position is uniformly zero, as appropriate for the 
disordered spin singlet, except for an abrupt dip centered 
at z = s Q . The dip is identified with a localized spinon 
at s. Therefore, ^ s s p may be thought of as the wave- 
function for a localized spinon at s. Starting from such 
an interpretation, in Ref.( M) we showed that, although 
spinons are collective excitations of a strongly correlated 
system, they can still be treated as real quantum me- 
chanical particles. In this paper we will generalize our 
formalism to states where both spinons and holons are 
present. 



The energy eigenvalue is 



^> v ] .:: = {-j(%i)( N -jf) 



B. The Norm 

The squared norm of the one-spinon energy eigenstates 
is defined as the scalar product: 
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E 



\^(zx,...,zm)\' 



(25) 



For the Haldane-Shastry model, we derived the formula 
for Eq.((25|) in Ref. §. By employing a recursion relation 
between (*^|^) and (*^_ 1 |* s ^_ 1 ), we expressed all 
the norms in terms of m and of the constant Cm intro- 
duced in Eq.fQ). The induction relation is: 



\ ^ r 



m(M 



(26) 



That recursively gives: 

T[M + l]T[m- 



/\& s p|\I/ s p\ 

\ m I ml 



T[M 



r[i]r[M + i]r[m + i]r[M 



(27) 



number but a coordinate variable. Therefore, unlike lo- 
calized one-spinon eigenstates, v f , ^° takes a well-defined 
crystal momentum, as we will show later. 

^>n° is a spin-singlet state. Indeed, by definition its 
total component of the spin along z is zero. Following 
the same steps leading to Eq.(p^), we also get: 



S z ^° = S 



= o 



which proves that is a spin singlet. 



h 



(29) 
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V. ONE-HOLON WAVEFUNCTION. 
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Holons are charged, spin-0 elementary excitations of 
the Kuramoto-Yokoyama Hamiltonian. They are con- 
structed by removing an electron from the center of a 


<S > 

ho 

-0.5 
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spinon. The state for a localized holon at ho is given by 
Ch i\^% ho ), where is the state defined in Eq.©. 



By construction, in the KY-model, the holon is the su- 
persymmetric partner of the spinon. However, unlike in 
the spinon case, the Brillouin zone for one-holon states is 
not halved, as both negative and positive energy holons 
can be constructed. In this section, we are concerned 
mainly with holon "kinematics" . We derive the one-holon 
eigenstates, their norm, their energy and their crystal 
momentum. In particular, we focus on negative-energy 
one-holon states, since these states are the ones relevant 
to the spinon-holon interaction. In the next Sections we 
analyze holon and spinon dynamics - the interaction be- 
tween spinons and holons and its relation to the instabil- 
ity of the hole excitations in the KY-model. 



A. One-Holon Spin Singlet 

Let N be odd and M = (N — l)/2. The wavefunction 
for a propagating, negative energy, holon is given by: 

M M 

^°( Zl , ZM \h) = (hr nfe - h >j life ^ zfe ) 2 > 

(28) 

where z\ , . . . , zm denote the positions of the | sites and 
h denotes the position of the empty site, all others being 
|. Also, < n < (N + l)/2. Different from the spinon 
case, in Eq.(E8h the holon coordinate h is not a quantum 
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FIG. 1. Spin and charge profiles of the localized holon 
defined by Eq.@. 6 is defined as 6 = -iln(z/h ), 
where z is the independent variable. 



B. Negative One-Holon Energy Eigenstates 

is an eigenstate of Hky with energy eigenvalue 
given by: 

WKY|^°) = {-J(^)(iV-l) 



J ,2-k 2 . N + 1, , . ,. 



(30) 



where < n < (N+l)/2. 

In order to prove Eq.(30), let us first split Hky as 
follows: 



J^ = h T s + h v s 

2 V N I 



IN 



l Q ^ Q 



(31) 



We define the various terms when calculating their con- 
tributions to the total energy. 

• Spin-exchange term: 

[h^°}(z 1 ,...,z M \h) 



6 



a#/3 



1 2a ™ 



1-iV 2 
24 



,1/ 



32 ~| C n.hor 



- A »+2 ^ U tt°°Czi,...,ZM|/Q 
9ft 2 J 1 ft" 

• Spin-Ising term: 

[ft£^°] (*!,..., 2 M |ft) 



(32) 



M M 

U. _ 7 .12 U, _ 



i=^3 J 



\ Zj h\* 



iV ( iV2 - 1 )\ 1 T,ho 



48 



^{ Zl ,...,z M \h) . (33) 



• Electrostatic repulsion energy: 

[h N ^l°}(z 1 ,...,z M \h) 

1 



-P 



(n Q + Ti/3) 



1 3 

4^-4 



P*^](zi,...,«m|/i) 



■("-charge kinetic energy: 

4 



[hU h n °}( Zl ,... : z M \h) 



E 

L 



T C /3T 



ho 



(zi, . . . ,2m|M ■ (36) 



To properly work out the contribution in Eq.(|36j), 
we have to express in terms of the .[-spin coor- 
dinates, 7/1, ... , rjM- In Appendix A we prove that: 

*£° (*!,..., z M I ft) = (%,... , 77m I ft) • (37) 

Therefore, we obtain: 

[h} Q ^]{ Zl ,...,z M \h) 





~N- 1 


h n - 






[ 2 "J 



A 2 - 1 77(77 - A) 
12 + 2 



1 +2 a 2 U j^^ 

2 a/7 2 / \ ft™ 



By using the identity: 



9ft 



-,n+l ^ 

9ft 



(38) 



A-l_ v _ft liT(]i „ 
9 ^ ft - 



<°(zi,...,Z M |/7) 



1 — A 2 
24 



V«°( Zl ,...,z M \h) . (34) 



• |-charge kinetic energy 

i , T ,ho 



[h^»°](z 1 ,...,z M \h) 



E 

*±P 



PCc 



\z a - zpl 



ho 
2 1 n 



(zi, . . . , 2Af|ft) 



v ^z^zp-hf d_ k ^l°( Zl ,...,z M \h) 
^ ^ k\\zp-h\ 2 W 

z jth k=0 1 p 1 

N 2 - 1 77(77 - A) 



h n - 


' N — 1 


77 




2 
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V+2 d 2 U ^°(z 1 ,...,z M |fe) 
2 9ft 2 Jl ft" 



9ft 



(35) 



(39) 



and the identity: 

un+2 d f^°{rn,...,T lM \h) 
9ft 2 



M 



(N-l)(N-2) | 2 g,/ ft 



ft — z. 



ft" 



we finally derive: 



ft™ 



(40) 



7 



[hl j ^ )(z 1 ,...,z M \h) 



M 



Oh 



~l hn '' 2 ~&\ v : " !/ h] • 14:1 ! 



Adding Eqs. (|32jj3^j3J,g5|,^ip all together, we find that: 
[H KY ^ }(z u ...,z M \h) 



J/27ry/ N(N 2 -l) 



2 \N 



48 



(42) 



+ n{n-M-l)j^°{ Zl ,...,z M \h) 
This is the formula we had to prove. 

C. Crystal Momentum 



The state \9„) is a propagating holon with crystal 
momentum 



7T , T 27T, 

-N H (n 

2 N y 



1. 



(mod 2ir) , 



(43) 



with the definition 
^°(ziz, . . . , z M z\hz) = exp(i 9 £°)^°(zi, . . . , z M \h) 



(44) 



Rewriting the eigenvalue as 



we obtain the dispersion relation 



I* 



ho\ 



(45) 



E{C) = -~ 



(f) 2 -(^°) 2 



(mod 7r) , (46) 



We worked out Eq.(46) for the case of negative-energy 
one-holon eigenstates. Unlike in the spinon case, the 
holon dispersion relation around the band minimum is 
quadratic in q, while it is linear in q near tt/2, where the 
band closes. 

Positive-energy one-holon eigenstates may be con- 
structed by supersimmetrically rotating one-spinon 
eigenstates Jig ]. Their energy is given by |_E(qJj°)| and 
the momentum spans the remaining half of the Brillouin 



zone. The corresponding dispersion relation is plotted in 
Fig. ([!]). Since for the purpose of studying spinon-holon 
interaction we only need negative-energy holon states 
only, we do not discuss here positive-energy holon states. 
We also need the wavefunction for a localized holon at 
site ho, which is obtained by Fourier transforming 

the propagating holon wavefunction to real space: 



(JV+l)/2 

E 

n=0 



. -ri^ho 



(47) 




71 


-7i/2 


7C/2 7C 




q 


N=9 








\-m— 


• — • 1 


— •! — • — *H 






N=ll 




•] — • — • — * — • 


• — !-• 1 



-7t -K/2 71/2 K 

FIG. 2. Upper graph: Holon dispersion relation. Positive 
energy holons are unstable towards decay into spinons; Lower: 
Configuration space for N=9 and N=ll sites 



D. The Norm 

The norm of the state 4'^° is defined as: 

<^°l*n°> = E \*n(Zl,...,Z M \h)\ 2 (48) 

From the definition of the norm it immediately follows 
out that (^"l^n } does not depend on n. The formula 
for (\l/^ |\l/^ ) is worked out in Appendix B and is given 
by: 



/*ho|*hov = ArM+1 (2M)!r[|]r[M + i] 

™ 1 2 M T[M + i] 



(49) 



8 



VI. ONE-SPINON ONE-HOLON 
WAVEFUNCTION 

As for many-spinon configurations, spinons and holons 
maintain their integrity when many of them are present 
in the same state. Therefore, we can diagonalize Hky 
within subspaces with a fixed number of spinons and 
holons. 

In this section we derive the action of Hky on the 
one-spinon one-holon eigenstates, diagonalize the corre- 
sponding matrix and work out the norm of the states. 



A. Action of Hky on one-spinon one-holon states 



• The Coulomb potential term acting on the one- 
spinon-one-holon wavefunction reads: 

[h N ^](z u ...,z M \h) = ^-^ . (53) 



• The | spin contribution to the charge kinetic energy 
gives: 

[h l Q ^](z u ...,z M \h) = 
- E ( ■■-.*m|z,9) 



To construct one-spinon one-holon eigenstates of Hky 
we start with states in a mixed representation, , where 
the spinon is localized at site s, but the holon is propagat- 
ing with momentum q^°. For N even and M — N/2 — 1, 
we have 

M M 

. . . , z M \h) = h n JJOj - s)(zj - h)zj Y[(zj - z k ) 2 , 
j j<k 

(50) 

where 1 < n < M + 2. To derive the action of the KY- 
Hamiltonian on , we split it as in Eq. (|3l"|) . 

• The action of spin-exchange term on \&™ provides: 

[h%*?](z u ...,z M \h) = 



EE 

3 = 1 zp^Zj 



2 

Z 3 Z 



;^{zi,...,Zp,...\h) 



N 2 -l n{n-N) 



N - 1 



12 



\Er« 

8 

dh V h" 



2 dh 2 \h n 
(54) 



• In order to manipulate the t~ s P m contribution 
to the charge kinetic energy, we need to express 
^(zi, . . . , Zu\h) in terms of the | spin coordinates 
rj. As in the one-holon case, we obtain: 

^( Zl ,...,z M \h) = % l {r ]l ,...,i lM \h) . 

Therefore, we have: 

[h^:]( m ,...,r l M\h) = 



+(M- 1) 



ds s + dh U" 



ds 2 s dh 2 \ h" 



1 / /i + s 
+ 2 I /i- s 



9s s 9/i U" 



(51) 



E 



[4^]( Zi ,...,^m|/i) - 

1 M N 2 -1 ^ 1 



12 



(52) 



{h-ztf 



S-/l) 2 

7V 2 -1 n(n-N) 



12 





~N — 1 
n 






ives: 


2 







2 d/i 2 I h n 



sh 



(s h} 



:^(Vi,---,Vm\s) 



(55) 
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By going back to f-spin coordinates z, Eq.(|55|) pro- 
vides: 



n 2 - 2n + 1 



1 — n \ h + s 
h — s 



hs 



2 s 



2 V h- s 



dh V ft n 



dh 2 V /i n 



n-l 



V|«,--fr| a * (.s-fr) 2 U/ h 1 j 



Adding up all the contributions together, we obtain: 
h KY ^(z 1 ,...,z M \h) = 



-iV 3 + 19JV / N , , , T 

h n 1 C 

48 V 2 



9 2 



(M- 1)3 — *™ - S^TTT^ 



9s 



9s 2 



\h + s f <9 T „ . T . 



/v,s 



(^- S ) 2 



*"-(s) 



n-l 



(57) 



-TV 3 + 19JV 
48 



On the other hand, the diagonalization of the "interac- 
tion" term, given by the fourth and the fifth rows of 
Eq. ( |57| ) , needs further work. By using a straightforward, 
although tedious, application of basic identities proved 
in the Appendix of Ref. [9j, one obtains: 



1 fh+s 



^ \1\h-a 

ses N 



+ (1 - nW s 



(h-s) 2 



1 h 



s^h 



1 fh + s 



(h-s) 2 



2 \h- 



lim s 

s — >h 



1 fh + s 



2 V h 



Sg~ s ^ + (1 - n)*? 



s—*« + (1 - n)** 



(h-s) 2 



(m - ?i + l)* mn + ^ iV(m - ra + l)* m _ i>n _j 
j'=o 



(60) 



In the next subsection we will will solve Eq. (p7j) by work- 
ing out the basis of one-spinon one-holon states that di- 
agonalize Hky- 



B. One-spinon one-holon energy eigenstates. 



if ; 



1 > and 



N 
~2 



M- 



(n - to - l)* m „ + ^2 Nin-m-l^m+j^+j 



(61) 



To diagonalize hx y, we introduce the propagating one- 
spinon one-holon energy eigenstates 

^ mn (z u ...,z M \h)^^2—^ n s {z 1 ,... 1 z M \h) . (58) 



if n — to — 1 < 0. Therefore the action of /iky on $ m „ is 
given by: 



N 

eo + rn(M — to) + n(n — 1 — ) 



The second and the third rows of Eq. (J57|) are diagonal 
in the basis of the states ty m n, and their contribution is 
given by: 



N 

[e + n{n - 1 - — ) + m(M - m)]* r 



(59) 



where: 



^ m 

-(n - to - l)*„ m - (n - to - 1) ^2 *m-j,n-i, (62) 



if to — rt + 1 < 0, and 



eo + m(M — to) + 7i(n — 1 — — 
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M—m 

-(m-n + l)* mn - (m-n + 1) ^ * 



(63) 



ifm — n + 1 > 0. Since the Hamiltonian matrix is upper 
(or lower) triangular, complete diagonalization is possi- 
ble. The energy eigenstates will be linear combinations 
of <3> m „ 



M—m 
3=0 



if to — n + 1 > 0. The coefficients are given by: 



6, = 



r[|]r[j + i] 



(69) 



(70) 



3=0 

if TO — H + 1 < 0, and: 

M—m 

*^mn — ^ '^m+j.n+j : 

J=0 



(64) 



(65) 



if to — n+l > 0, corresponding to the energy eigenvalues 



m " 2 V N 



eo + m(M — to) 



TV 1 

+n(n — 1 — ) — — (n — m — 1) 



if to — n + 1 < 0, and: 



e: 



J /2tt 



2 \N 



€q + m(M — to) 



AT. 1, 
+n(n - 1 - —) + -{n - m- 1) 



if m — n + 1 > 0. The coefficients a/ are defined by the 
recursion relation: 



ai 



1 

2/ ^ 

k=0 



a k a Q = 1 



(66) 



In terms of spinon and holon momenta, E^ nn and E mn 
take the same form E mn , given by: 

Emn = E GS + E(q%) + - ^ |g " = . (67) 

i? mn is the sum of the ground-state energy, the energies of 
an isolated spinon and an isolated holon plus a negative 
interaction contribution that becomes negligibly small in 
the thermodynamic limit. Eqs.(|64|, |65| ) can be inverted. 
The result is: 



3=0 



(68) 



if to — n + 1 < 0, and: 



C. The Norm 



The squared norm of the state ^> m n is defined as: 



Z1,...,ZM 



|$ mn (zi, . . . , z M \ 



(71) 



In a similar fashion to the two-spinon case discussed in 
we compute the norm of the one-spinon one- holon 
states by means of mathematical induction. The calcu- 
lation is presented in detail in Appendix C. The basic 
induction relation in the case n — A; + 1 < is given by: 



kn \ ^kn/ 



<*fe— a,7*|*fe— x,- 



(fc- l)(M-fc + f) 
k(M -k + 1) 



(72) 



which provides the formula for the norm of the one- 
spinon one-holon energy eigenstates: 



($ kn \<S> kn )=N M +'-^f(M+^ 



. 1 (2M)I 



i.r[M-* + i]r[* + i] 



2 y r[M-fc + f]r[fc + i]' 

(73) 

In the complementary case, n — k + 1 > 0, we obtain: 

(®kn\&kn) = (®M-k,M-n\®M-k,M-n) 



i r[M-fc + |]r[fc + i] 

2^T[Af - A; + l]r[fc + f] 



(74) 




FIG. 3. Left panel: \p mn {e lB )\ 2 vs. 9 for m = M,n = 0; 
Right panel: the same plot on a log-log scale. The dashed 
straight line is a plot of 1/9 
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VII. SPINON-HOLON ATTRACTION 

In this Section we analyze the interaction between a 
spinon and a holon by constructing the real-space repre- 
sentation of the one-spinon one-holon wavefunction, and 
by studying the behavior of the corresponding probability 
as a function of the separation between the two particles. 
Our exact results show that a spinon and a holon interact 
through a short-range attraction identical, in the thermo- 
dynamic limit, to the atraction between two spinons. 

The state for a localized spinon at site s and a localized 
holon at site ho, ^ s h is defined as the Fourier transform 
of <3>™ back to coordinate space: 



M+2 



= E v*i 



(75) 



Following the same steps as for the two-spinon wavefunc- 
tions we define the real-space coordinate represen- 

tation for a spinon-holon pair, s m h ~ n p mn (s/ho), as fol- 
lows: 

M+2 n-2 m 

71 — 1 771 — j—0 



M+2 M M-m 
E E S " l V l E b ^rn +j ,n +j 

n—1 771— 71 — 1 J—0 
M+2 n-2 



= E E * m vw^)<£ 



n— 1 m— 



ho + S 

h - s 



h 



dhn 



shot 



- ho 



(78) 



where v = M \im — n+1 < 0, v = otherwise. In the 
differential operator in Eq.(|7^), we recognize the sum 
of the energies of the free spinon and holon, a velocity- 
dependent interaction, which diverges at small spinon- 
holon separation, and another term which takes into ac- 
count the correction for the case when the spinon and the 
holon are at the same position. By using Eqs. (|7^ , |77| . |78| ) . 
we find the following equations for the "relative wave- 
functions" p m n(z) and p mn (z) (z = s/ho): 



[2- 



1 



dz (1 — z) 
if m — n + 1 < 0, and 



]Pmn(z) + 



y M-m-l 



Pmn( Z ) 



(1-z) 



(1-1) 



-2W(1) = , (79) 



p mn (l)=0, (80) 



if m — n + 1 > 0. 

Eqs.( f79| , [j0"| ) are first-order "Dirac-like" equations. 
They are first order because the spinon and the holon 
energy bands have opposite curvature. In this respect, 
they differ from the differential equation obtained in the 
two spinon case, which was second order (8|J^]. The cor- 
responding solutions are given by: 



Pmn {%) 



M-m-1 

E 

fe=0 



r[fc- 



r[i]r[fc + i]' 



(81) 



M+2 M 

E E S m ho n Pmn (^r, 
n=l m—n—1 



1*0 



(76) 



where | $ m „) is an eigenstate of Hky with eigenvalue 

J mn, 



($nrn|#KY|* s /w) = E mn {$ mn \V sho ) 



(77) 



The matrix element {<&mn\HKY\^ sh ) can be written as 
a differential operator acting on the analytic extension 
of (Qmnl^sha), where s and ho are understood to take 
any value on the unit circle. Therefore, by equating 
{^m-nXH-KY^ sh Q ) to E mn ($ mn \$ sho ) , it is straightfor- 
ward to write down the equation of motion for the one- 
spinon one-holon wavefunction, which reads: 



[E mn - E GS )(<5> mn \^ s 



($ m n\(HKY-E GS )\V sho ) 



d d d N d 

{M ~ s lTs>irs +h °Wo {1 + Y + h °dh-o ) 



for Eq.fl79j) and 

Pmni z ) 



E 

fc=0 



r[fe- 



r[i]r[fc + i] v z 



(82) 



for Eq.(gO|). 

The value of the spinon-holon wavefunction at zero sep- 
aration between the particles is derived in Appendix D. 
It is given by: 



and 



Pmn(l) 



Pmn ( 1 ) 



r[M- 



T[l}T[M-m] 



'r[i]r[m + i] 



(83) 



(84) 



Within the framework of our formalism, it is possi- 
ble to treat spinons and holons, collective excitations 
of strongly-correlated one-dimensional electron systems, 
as actual quantum-mechanical particles. We were first 
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able to associate a two-particle wavefunction to a spinon- 
holon pair and to write down the corresponding equation 
of motion (Eqs.(76 7S|)). We then worked out the exact 
wavefunctions corresponding to each energy eigenvalue 
(Eqs.(|l]fl)). 

The squared modulus for the spinon-holon wavefunc- 
tion, \p rnn (z)\ 2 , gives the probability for a spinon and a 
holon configuration as a function of the separation be- 
tween the two particles. In FigJ| we plot |p„ m (e* e )| 2 
versus the distance between the spinon and the holon, 9. 
From Fig|| the nature of the interaction between a spinon 
and a holon may be easily inferred. While at large sepa- 
rations \p m n{e ie )\ 2 does not depend on 9, as it is appro- 
priate for noninteracting particles, at small separations 
it shows a remarkable enhancement. This corresponds to 
a huge increase of the probability of configurations with 
the spinon and the holon on top of each other. The en- 
hancement does not depend on the holon momentum. As 
N gets larger, the probability enhancement peaks up. It 
survives the thermodynamic limit, even though the inter- 
action energy goes to zero, and the total energy becomes 
the sum of the energies of the isolated spinon and holon. 
However, the attraction is not strong enough to create 
a spinon-holon bound state, even in the thermodynamic 
limit. This corresponds to the absence of a low-energy 
stable hole excitation, and it is what causes the quasi- 
particle peak to disappear. 

An intriguing feature of the spinon-holon interaction in 
the thermodynamic limit is that it has the same power- 
law form as the spinon-spinon interaction derived in [^|,^| . 
In the right panel of Fig.[| we plot \p mn (e l8 ) | 2 on a log- log 
scale and compare it with 1/9. The probability falls off 
as the first power of the separation between the two par- 
ticles. This shows that, although the equation of motion 
for a spinon and a holon is quite different from the one 
for two spinons, the interaction in both cases results in a 
short-range attraction, and its effects on the correspond- 
ing two-particle wavefunction are basically the same. 



VIII. HOLE SPECTRAL FUNCTION 

In this Section we work out ^ P:ho (a;, q), the one- 
spinon one-holon contribution to the hole spectral func- 
tion Ah(tu,q). We show that this contribution (which 
provides quite a good approximation to Ah_(to,q) for 
q ~ 0) depends only on the p m „'s and the p mn 's cal- 
culated at z = 1 (that is, as the spinon and the holon lay 
at the same site). This allows us to obtain for any finite 
N a simple closed-form expression for A^ p ' ho (a;, q), and to 
relate it to the spinon-holon interaction. In the thermo- 
dynamic limit, we obtain the previously known formula 
for the contribution of the one-spinon one-holon states 
to A^' °{u>,q) 14|. The formula in |u| shows that there 
is no low-energy hole pole in the hole spectral function, 
but rather a sharp square-root singularity followed by a 



branch cut. These features have also been experimen- 
tally detcted by means of ARPES experiments on quasi 
ID insulator @. 

The branch cut corresponds to the lack of integrity 
of the hole excitation, which breaks up into a spinon 
and a holon. Here we will show that, in the ther- 
modymanic limit, the probability enhancement p mrl (l) 
(p„ m (l)) turns into the square-root singularity at thresh- 
old for a spinon-holon pair. As a consequence, we prove 
that the square-root singularity in the hole spectral func- 
tion is a direct consequence of the interaction between 
spinons and holons. Therefore, it can be directly experi- 
mentally measured. 

We begin with the calculation of A^ p,ho ((jj, q) for a finite 
lattice. In Lehman representation we obtain: 



A^' ho (u), q) — 3m 



E 

X 



l(*IE, l0 (M- fc Q i0 Tl*Gs)| 



TN(X\X)(y GS \V 



GS) 



1 



uj + - (E x - Eqs) 



(85) 



where \X) is an exact one-spinon one-holon eigenstate of 
Hky, \®mn) w hh energy E x = E mn , and q = 2nk/N 
(below we discuss why only forward propagating states 
contribute to Eq.(|S5|)). 

Using v4^ p ' ho (w, q) instead of Ah((v,q) is equivalent to 
approximating 



GS/ 



hoho 



M+2 n-2 

V V ft?- n p mn (i)$ B 



n— 1 rn—0 



M+2 M 

E 

7i=l m=n—l 



E K-y'mnim. 



(86) 



(Eq. (|6|) basically amounts to neglecting contributions to 
c/iotl^cs) coming from multi-spinon one-holon states.) 

Since Hky contains the Gutzwiller projector P, its ma- 
trix elements between states with at least a doubly oc- 
cupied site are zero. Therefore, at half-filling, Ah(q, u) 
takes contributions only from forward-propagating hole 
states. Hence, using Eqs.(|85|,|8^) we obtain: 



,»sp,ho/ \ 



1 (M+2 1-2 
k 1=2 m=0 



6 k -m+ip 2 m l(l) 

■ ir) - (E mn - E GS ) 



M+2 M 

T T 

1=1 m=l-l 1 y 



S k - m +l(p ml ) 2 {l) \ ($ m l\$ml) 



E G S) J (*GS|*GS> 



(87) 



Eq.(p7[) shows that only the p mn 's at z = 1 determine the 
spinon-holon contribution to the hole spectral function. 
Therefore, the contribution is completely determined by 
the spinon-holon interaction. 
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Let us now analyze the thermodynamic limit of 
Eg. (|S7|) . In the thermodynamic limit, the gamma func- 
tions can be approximated by using Stirling's formula 

T[z] pa V^(z - l) z -h-( z -V 

From Eqs.( |S7| , |S8| ) we get, in the thermodynamic limit 

* sp.ho 



^' fa") 



7f(M + 1) 



■ M+2 1-2 



r— m+/ 



'Af -m - i 



E G s) 



rr> 



M+2 M 

E E 

Z=2 m=/-l 



(5 r - 



■m-\-l 



lu + irj - {E mn - E GS ) V M - m 



(89) 



(Notice that, in order to stabilize the hole occupation at 
1, we had to introduce a chemical potential Jit 2 /A, which 
is added to the energies E mn .) 

By defining the auxiliary variables: 

2tt 2tt 

Eq. (|S9"|) , in the thermodynamic limit, may be written in 
the form already obtained in |Q: 



9ho dg sp / 7T - g sp 

7T V 9sp 



IX. CONCLUSIONS 

In this paper, we have extended the formalism in- 
troduced in H to analyze spinon interaction in the 
Haldane-Shastry model, to the case where also charge 
degrees of freedom are involved. Our formalism allows 
us to define a quantum-mechanical real-space represen- 
tation of the one-spinon one-holon wavefunction. Wc 
construct a Dirac-like equation, whose solution is the 
spinon-holon wavefunction in real space coordinates. By 
means of a careful study of the real-space one spinon one 
holon wavefunction, wc show the existence of the spinon- 
holon interaction and its survival in the thermodynamic 
limit. Spinon-holon interaction generates a short-range 
enhancement in the probability for a spinon and a holon 
to be on the same site. The attraction, however, is 
not strong enough to form a spinon-holon bound state, 
which would correspond to a Landau's quasihole reso- 
nance. This makes, in the thermodynamic limit, the 
hole excitation fully unstable against decay in one-holon 
multi-spinon states and the quasiparticle peak disappear. 
Correspondingly, in the thermodynamic limit, the prob- 
ability enhancement develops into a square root singu- 
larity followed by a branch cut, which reflects the full 
instability of the hole excitation. Hence, by means of a 
sequence of exact, straightforward steps, we prove that 
spinon-holon attraction is what makes Landau's Fermi 
liquid theory break down in 1-d strongly correlated elec- 
tron systems. 
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dq sp / gs P \ S(q-q sp + q ho ) 
,,.,„ 7T y TT-q sp ) uj- fi + irj- E(q sp ,q ho ) 



(90) 



In the region < q < n the integration of Eq.(|90|) gives: 



K Ph °(^q) 



1 J[g + f] 2 



n 2 JqV CJ-Jfe-f] 2 
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J[— + q(^ -q)]-u 



(91) 



This formula shows that the spinon-holon probability en- 
hancement p^ n (l) turned into a square-root singularity 
in A^' °(ui,q) at the threshold energy for creation of a 
spinon-holon pair. Because the spinon-holon joint den- 
sity of states is uniform, the main conclusion we trace 
from our calculation is that the sharp nonanalytic thresh- 
old in A^' °(ui, q) is the direct consequence of spinon- 
holon interaction. 



APPENDIX A: FUNCTIONS EXPRESSED IN 
TERMS OF j-SPIN COORDINATES. 

In this Section we prove the formulas that express the 
states of the KY-model in terms of the J,-spin coordinates, 
once the expression in terms of the | _ spin coordinates 
is known. The starting point is the following identity 
(z Q , zp are N — th roots of the identity): 



n ^ 



lim 



1 



♦2/3 Z— Zfi 



N 

Z/3 



(Al) 
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The ground state wavefunction at filling- 1/2 expressed in 
terms of the "f-spin coordinates is given by: 



At 



A I 



^ G s(z 1 ,...,z M )=Y[(zi~z J ) 2 Y[z t , (A2) 



where N is even and M = N/2. Let 771 , . . . , r\M be the 
J,-spin coordinates. Upon applying Eq.(Al), we get: 



MM MM 

n> - z 3 f n zt = (-1)"^ - ^) n ** 

%<j t i=£j t 



(-1) 



MiM+nllt Zt iit \j:> 



M 



i<3 



(A3) 



which proves that ®gs(zi, . . . , z M ) = ^gs(Vi, Vm) 
The one-holon wavefunction is given by: 



M 



M 



M 



. . . , ZM \h) = jjfe - ft) - n z * > 



(A4) 



where now N is odd, M = (N — l)/2 and ft is the coor- 
dinate of the empty site. 

The same steps as for *&gs a Pply to the one-holon 
wavefunction. We have: 



M 



M 



M 



A I 



A I 



M 



ft"(-l) 2 



^ ^"3 J \_\ z i 

i<3 I 



^_±_i> fit z t llf'djr) 



Y[ Vi , Zj {zj - m) 



AI 



A I 



A I 



h n H(Vi-s)(m-h)l[( Vi - Vj fl[ m . (A7) 



Eq.(A7) provides the proof that ty™(zi, . . . , Zm \h) = 

The last identity we need refers to the case where the 
spinon and holon are at the same site, which we had to 
consider in deriving Eq.mq). We have: 



A I 



A I 



M 



. . . , z M \s) = s n Y[( Zj sf Y[( Zl - Zj f [] zt 



M 

n 



z<3 



(z 3 - S ) gif^xnng)) 

(zj-h) 11. .,, (zj //;) 



fr fe-.s) (nf ^)(nf(f)) 
A /(*i-*) nf(f) 



M 



M 



i<j i 



ft"(-l)- 



. lit z *llt (-) 

n m ,zi(zj - vi) 



A I 



A I 



A I 



A I 



A I 



i<j 



AI 



At 



AI 



i<3 



(A5) 



this proves that *£°(>i, . . . , z M \h) = ^°(r?i, . . . , rj M \h). 

The one-spinon one-holon state ^"(^i, . . . , Zjvf|ft) is 
given by: 

*™(z 1 ,...,z M |ft) = 



h n H(z ] - S )(z ] -h)Y[(z l ~z ] fl[z t , (A6) 



where N is even, M = N/2 — 1, s is the coordinate of the 
J,-spin and ft is the location of the empty site. As in the 
previous cases, we have: 



= -(tY 1 - 1 ^ Him - hf Vj f Yl th . (A8) 



Eq. (|A8|) proves the identity 



^( Zi ,...,zm\s) = -(£)" 1 n(vu---,m\h) 



APPENDIX B: THE NORM OF ONE-HOLON 
WAVEFUNCTION. 

In this Appendix we discuss in detail the calcula- 
tion of the norm of the negative-energy one-holon states, 
(**°|**°). 

(*£°l*£°> = E l^°(Zl,---,*M|ft)| 2 • (Bl) 

Zl,...,ZM,h 

Let P m (zi, ■ ■ ■ , zm) be the m — th degree symmetric poly- 
nomial in zi, . . . , zm- One obtains M: 
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*™(zx,...,z M \h) 



A I 



M 



= E h m+n P m ( Zl , z M ) Y[(zi - zj) 2 



m=0 



From Eqs.(Bl,B2), we derive 



M 



m— zi,...,zm 



M 



(B2) 



of the norm of one-spinon one-holon energy eigenstates. 
As in the calculation in Refs. |^,^|, the key operator is 
given by ei(zi, . . . , zm), defined as: 



ei(zi, . . . , zm) = z\ + . . . + z M 



(CI) 



The state for one holon and one spinon localized at s is 
given by: 



z M ) = $"(21, ... , z M \h)^Gs(zi, . . . , z M ) 



(C2) 



,1/ 



(B3) 



m=0 



The norm of one-spinon wavefunctions, (^ml^m), has 
been derived in Ref. [01, where it has been shown that: 



r[M + i] r[m + ^]r[M-Ti 

r[M+ |] T[m + l]T[M-m + l] ^ 2 M 



(2M)\ 



2i N M\~^r^ (B4) 



Therefore, we can re- write Eq.(B3) as 



' 71 I 71 



r[ 



2 J 2 



M 



5mo(1) 



(B5) 



where the polynomial g m n(z) is the two-spinon relative 
wavefunction, denoted by p mn (z) in Ref. M: 



Qmn (z) 



T[m — n + 1] 

r[|]r[m - n + |] 



where: 



a/ 



$:(z 1 ,...,z M |/ l )=/ l "nfe- s )fe-M ■ 

On ^"(zi, . . . , zm), ei acts as a ladder operator, as we 
are going to show next. 

In order to work out the action of Hry on holon eigen- 
state, we splitted it into five terms: Hky/^{^-) 2 = 
tig + tig + /jjv + hq + hq. Among those five terms, the 
only ones that do not commute with e\ are the spin ex- 
change operator tig and the |-spin charge propagation 
operator tig. On the state 'J™, tig is realized as: 

f 1 M d 2 

h T s^ n s = * GS {const. + - zf-g^ 



M n 

zf 



d A — 3 



Z% Zj 3zi 



t * 



(C3) 



E 

fc=0 



r[fc + ^]r[m-n-fc + ^] fc 
f[jfe + l]T[m- n - jfe + 1] Z 



(B6) 



By means of generic properties of hypergeometric func- 
tions [pTgfl , one gets: 



Qmn ( 1 ) 

Therefore, we obtain 



r[^]r[m-n + 1] 

r[m — n + ^] 



(B7) 



( ^ 0|e) = iV M +1 (2M)!r[i]r[M + i] 
2 M r[M + §] 



From Eq. 



we derive that: 



GS\ 



M - 



er^+E 



* GS <j (A/ + -)e 1 ^ + M( S + / l )^ 



ds 



(C4) 



APPENDIX C: THE NORM OF ONE-SPINON 
ONE-HOLON ENERGY EIGENSTATES. 

In this section we generalize the recursion procedure 
introduced in Ref. J8||| to calculate the norm of one- 
spinon and two-spinon wavefunction to the calculation 



In the sense clarified by Eq.(C4), e\ commutes with tig , 
hff and tig. On the other hand, it does not commute with 

tig. Indeed, since, in order to derive the action of tig on 
we have to express the state in terms of the |-spin 
coordinates, we must do the same with e±. In order to do 
so, we notice that { Z j}, h and s taken all together 



1G 



are the set of the N n — th roots of 1. Therefore, we 
obtain: 

z\ + . . . + Zm + Vi + • • • + Vm + h + s = , 
which implies: 

ei(zi, . . . , z M ) = -ei(r?i, . . . , t)m) - h - s . (C5) 
The action of hg on ei\&™ gives: 



iV 2 - 1 n(n - N) 
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ei^(MIM 



iV-1 



(-/i-s-^?7j) 



ei ]*« = (M + i) ei *» + M( S + h)^ s 



J_ ( 2tt 

2U 



s 2 — + 



9s 



7 Af 

h \ - 

h-s ^ 

m=0 



(C8) 



where we have introduced the one-spinon one-holon plane 
waves ^m n defined in Section VI-B. 



In order to further manipulate term in Eq.(C8), let us 
consider, now, the identity: 



m=0 



1 B 2 

2 dh 2 



M 



(-*--£ 



*?(MI>0 



A/ 

E 

m=0 



h l+rn 



-k 



N 



E 



( S -/*) S 



E ^ ~ 2s 



n-l 



N 



*mnh 1+m - k E 



m=0 



n — m — 1 , , 

x - / s \ ™+r-fe 

{s/M r=0 



= *Gsei[/i^r](i7i,...,77M|/») 



JV- 1 



^(J7i,...,^|/i) 



3 f^( m ,...,r, M \h) 

dh 



M m—n—1 

i E E £ (|) 



n+r— fe 



(C9) 



m— n+1 



{s/h} r=0 



(The symbol X){ s //i} means that we have to sum over 

s /heS N ). 

Suppose, now, n — k + 1 < 0. In this case, the sum in 
Eq. (|C9|) over m from n + 1 to M will give 0, while the 
second sum will be reduced to: 



+ <t>^( m ,...,r, M \s) 

n — s 



Eq.(C6) yields the result: 



^gs< h 



(C6) 



\zi, . . .,z M \h)-y-j ^h(zi, . . .,z M \h) 



(C7) 



Upon summing Eq.@ and Eq.Q, we derive the ba- 
sic relation we need in order to work out the recursion 
relations: 



E" 1+m ~ fc E E 



m=0 



{s/h} r=0 



m+r — k 



On the other hand, as fc — n + 1 > 0, the only nonzero 
term in the sum will be given by: 



M 



E *mnh l+m - k E E 



m—k+1 



{s/h} r=0 



a \ n — k+r 



In order to write down the action of [i?KY,ei] on the 
pla ne- w ave states ^km we multiply both mem ber s of 
Eq. (|C7j) by s~ k and sum over s. By using Eq.(CE), it 
is straightforward to derive the following equations: 

• lffc-n+l<0: 

[*™,et]* kn = (M+he 1 * kn 

2^1 Z 
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+{M -k + l)* fe _i,„ + [M -n+ l)h9 



fc-i 
m=0 



• lffe-n+l>0: 



h l+rn-k ^ "^I^ 1 ^^m-fc+r 



s e a) 



{s/M r =o 



(CIO) 



2 V AT y 



ei]* fcn = (M + -)ei* fe „ 



n — m — 1 

N ^ ^ \h 

{s/h} r=0 



^ / e \ m — k+r 



(C12) 



From Eq.(pi3), we get: 



{E k -i, n - E kn -M- -)($fe-i,„|ei|$ fcn ) 



fc-i 

^a^(M-fc + £+l)($ A; — l,n 

*=0 



+(M - fe + l)**-i lf . + (Af - n)/i* fe „ 



h 1+ 



M 



^ ^ - ran 



m — k 



m —71—1 



{s//i} r =0 



(Cll) 



Let us now work out the basic recursion relation in 
both cases. 

• Case k — n + 1 < 0. 
In this case energy eigenstates are given by: 

k 



e=o 



(see Section VI-B for the definition of the coeffi- 
cients dj). 

Therefore, we have: 



Hky i*. r Hky i.t, 



2 V Ar J 



J/2tt\2 : 

£=o 2 v at ; 



k 1 

= Y j ae(M+-)ei$ k - 



1=0 



k-l 



^a*(M- k + i + 1)* 



fc-l-£,n- 



£=0 



■ ^ a*(M - n + 1 + l)^ k -t,n-i+i 



Yai(M-n+£+l){^k-i, n \^k-£,n+i-e)+{^k-i,n\^k-i,- 



e=o 



= [M-fc+2+6i(Af-n+l)+oi(Af-n+2)](*fc_i in |*fc_i, n ) 



= (Af - fc+ -)($fc-l,„|$fc-l. 



(C13) 



(notice that, from their definition, we have a\ 



-1/2, &a = 1/2). Eq.( |C13| ) may be recast in the 
following compact form: 



<* fc -i,„|ei|$ fcn ) _ M-fe + | 



(*fe-i,„|$k-i,») 2(M-A + 1) 



(C14) 



• Case fc — n + 1 > 0. 
In this case we have: 

M-k 



&kn — a£$fc+£ )fl +* 



t=0 



Eq.(C12) now takes the form 



J /2tt\2 
2 ^ AT I 



M-fe 



«=0 



M-k 



Y at{M -k-£+ l)y k+e - hn+ i 



e=o 



k-l k-t-1 



5Z ae ^m,n+l+m-Ai 



=0 m=0 



M-k 



,n+i+e 



e=o 
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M-k M 

^2 at ^2 *m ,n+l+m — k 
1=0 m=k+e+l 



APPENDIX D: SOLUTION OF THE EQUATION 
OF MOTION FOR THE ONE-SPINON 
ONE-HOLON WAVEFUNCTION 



N 



E 

{s/h} 



i— n— l 



n— fc+r 



r=0 



(C15) 



From Eq.(C15) and by working exactly as in the 
previous case, we get the identity: 

(E k -i, n - E kn - M- -)($fe-i,„|ei|$ fcn ) 

= (M-fc + l)($ fc _ 1 , n |$ fe _ 1 ,„) , (C16) 
which yields the relation: 

M - h -I- 1 

(C17) 



k-l,n|*k-l,n) 2(M-fc+i; 



In order to complete the recursion procedure, we 
need one more induction relation that is derived by 
inserting in the product (<fr a b|ei|3> c d) the product 
zm£*m = 1) where: 

fiM^ii ■ ■ ■ , %m) — z\ ■ . . . ■ z M ■ 



From the definition of scalar product, Eq.(|25|), it is 
straightforward to show that: 

($afc|ei|$ cd ) = ([(e M ) 2 $a b ]|ei|[(eM) 2 $ C( i]> 



(®M-c,M-d\ e l\®M-a,M-b 



(C18) 



By app lying Eq.( |C18| ) to Eq.( |C13| ) and by using 
Eq.( pl7 ), we obtain: 



($fc-l,n|ei|$fc„) = ($M-k,M-n\ei\$M-k+l,M-n) 

= - 2{k k _ 1 _ ) ^kn\^kn) ■ (C19) 

By putting together Eq. (|C15 ) and Eq.(C19) one finally 
obtains: 



($fcn|$fc™) 



(4> 



fe-l,n $k— l.n 



(fc_ i)( M -fc + |) 
fc(M- fc + 1) 



(C20) 



that provides the formula for the norm of the one-spinon 
one-holon energy eigenstates in the case k — n + 1 < 0: 



In this Appendix we will derive the solution to the 
equation of motion for the relative coordinate part of the 
spinon-holon wavefunctions, p m n{z), p mn (z). In order to 
do so, let us consider first the case m — n+1 < 0, in which 
we express the solution to Eq.(|79|) as a power series of z: 



Pmn(z) = 



a k z 



(Dl) 



From Eq.(Dl), we get the following equation for the co- 
efficients a k '. 



2^fca fe+1 z fe + ^(2fc-l) a k z 



w M~ra 



5> = o. 

k 

(D2) 



As k < M — m, the following recursion relation between 
the afc's holds: 



a k k + 1 
Eq. (|D3|) is satisfied by 



a k = 



T[k + ±] 

r[i]r[fc + i] 



(D3) 



(D4) 



To calculate au-m+i, we need the following identity, 
valid for any positive integer R (Co is a closed path cen- 
tered at z = 0): 



dz 1 



1 



f^ r[fc + |] A 

h r I & + !] h Jc a 2m z^ 



dz z R+1 - 1 1 
Co 2niz R+1 z- 1 y/1 - z 



dz 1 
Oo z R + l 2m (1 - z)f ~~r[i?+l] v " ' 2' 



T[R+h], 1, 



(D5) 



($kn\<S>kn)=N M +'±=^(M+-) 



-i(2M)! 



1 r[M — fc + l]r[fc + -] Since the recursion relation for aja-m is: 



2 m v— ■ 2 'r[A/ -fc + §]r[fc + i] 

(C21) 

In the complementary case, k — n + 1 > 0, we may follow 
the same step to prove that: 

($fcrJ$fc«) = (®M-k.M- n \®M-k.M- n ) 



Af , , (2M)' 1 . T[M — fc + h]T[k + 1] , 
= N M+1 ^-rJ-(M + -)-^— ' ^ ; 3 ■ (C22) 



-2(M - m)a M - m + 2(M - m - -)a M - ro -i 



(D6) 



Eq. (|D5|) implies ajvf_ m = and: 

M-m-l 



2 M 



2 y r[M-/c + i]r[fc + §] 



r[&- 



fc=0 



r[i]r[fc + i] 



(D7) 
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Eq.(D7) provides the formula for p mn {z) used throughout 
the paper. 

The same steps followed in the case m — n + 1 < 
allow us to find the spinon-holon wavefunction in the case 
m — n + 1 > 0, provided the coordinate z is substituted 
with £ = l/z. 
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